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SUMMARY 


This report contains the algorithms necessary for constructing an 
aircraft flight path which contains a glideslope change at a waypoint which 
lies along a straight line. The report also contains the necessary algorithms 
to reconstruct the glideslope change waypoint along a straight line in the 
event the aircraft encounters a valid MLS update and transition in the ter- 
minal approach area. 

Results of a simulation of the Langley B737 aircraft utilizing these 
algorithms are presented. The method is shown to reconstruct the neces- 
sary flight path during MLS transition resulting in zero cross track error, 
zero track angle error, and zero altitude error, thus requiring minimal 
aircraft response. 



INTRODUCTION 


Previous studies have shown the advantages of eliminating cross track, 
track angle and delta altitude errors through the method of a small path re- 
construction at transition from TACAN to MLS navaids in the terminal area 
(Refs. 1, 2, and 3). These studies have been confined to a path which 
allowed vertical glideslope changes to occur only at the centers of turns. 

Since some terminal area flight paths require glideslope change to occur along 
a great circle (straight line) path, this study was undertaken to include glide- 
slope change waypoints defined along straight line segments and to develop 
algorithms for updating these waypoints, at transition, in a manner consistent 
with the Continued Track Method outlined in Reference 3, in which the cross 
track error, the track angle error, and the vertical altitude error are all 
eliminated at the point of transition. 

The algorithms are written in vector notation, consistent with the 
guidance and RNAV equation used in References 1, 2, and 3. 

The report contains the derivation of the algorithms and the results of 
a computer study utilizing the equations for a variety of aircraft and runway 
landing conditions. 
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I. The Glideslope Change Waypoint Along a Straight Line 


The input data required to define a great circle straight line on the Earth’s 
surface consists of the latitude and longitude of two consecutive waypoints on 
that line. Thus, given Xj t ^I+l ’ we have 


and 


WR (1) = 



(6j) 

(6j) sin (Xj) 
(6j) cos (Xj) 


WR (2) 


sin 

<-cos(Vl)^‘n Vl> ? 

1 con nos 


( 1 ) 


( 2 ) 


A 

The unit normal, WN(I) , perpendicular to the plane of the great circle 

A A 

which contains WR(I) and WR(I+1) is given by 


WN(I) = 


WR(I) X WR(I+1) 

A A 

|wR(I) X WR(I+1)| 


(3) 


The difficulty in defining the input data for three consecutive waypoints, each 
of which is to be required to lie in the same great circle plane , is that one has 
to solve a nonlinear equation for the intermediate latitude and longitude of the 
middle waypoint. In order to simplify the input data requirements, it is 
sufficient to input an approximate latitude and an approximate longitude for 
the interior point not necessarily on the line, plus a numerical index identi- 
fying a glideslope change wa 3 q>oint. The path generation algorithm, stored in 
the computer memory, will automatically relocate the approximate wa 5 T>oint, 

A 

WR(I+1), associated with proper index, on to the great circle, utilizing the 

A 

cross track error algorithm for finding the desired unit vector, WD , on a 
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great circle path, given the unit waypoint, WR(I+1), and the unit normal, 

A A A 

WN(I), defined by the vector cross product of WR(I) and WR(I+2). 

The index recommended for use in this report is to define the radius 
of desired circle, R_(I), associated with glideslope change waypoint, 

A JL 

WR(I+1), to be zero. 

Thus, in computing the unit normal. Equation (3), between two successive 
waypoints, if the radius of the turn circle, R„(I), associated with the way- 

A J. 

point WR(I+1) is equal to zero (0. ), the algorithm for determining the unit 
normal is given by. 

If, R^(I) = 0. (4a) 


then 


W(I) = WR(I) X WR(If2) 


WR(I) X WR(I+2)| 


(4b) 


and the relocated glideslope change waypoint, WR(I+1) becomes 




WD = ^ ( WR(I+1) - sin b WN(I) 




J 


(5a) 


where 


cos b = WR(I+1) • WN(I) 


and 


(5b) 


sin b = \l 1 - ^ 


cos b 


Finally, if the radius of the turn, R (I) , associated with the glide- 

A X 

slope change waypoint, WR(I+1), is zero, we determine the unit normal, 

A A 

WN(I+1), to be equal to the previous unit normal, WN(I). Thus, 


Then 


WN(I+1) = WN(I) (6b) 

The only additional change that is required in the algorithm for path 


construction is to set the required turn angle, 0 (I) , associated with the 

A 

glideslope change waypoint, WR(I), whose radius, R^(I) , is equal to zero. 


Thus, 



If, 

R^m = 0. 

(7a) 

Then, 

o 

II 

(7b) 


This is required because the computation for the turn angle, 0 (I), is 
obtained by means of an arc tangent routine, and some round off inaccuracies 
could result in a setting 0(1) equal to some small angle, thus producing an 
error in the distance to go logic. 

The geometric equivalent of Equations (4b), (5a), and (6b) are shown in 
Figure 1. 

For a normal first interior waypoint, containing a finite non-zero turn 
radius, R^(l), the path generation algorithms will produce two contiguous 
segments. The first segment is a great circle (straight line) beginning at the 

A 

initial unit waypoint, WR(1), and terminating at the first incoming tmit tan- 

A 

gent vector, PI(1), which marks the start of the turn. The second segment 

A 

is the first half of the turn about the first unit center turn vector, CR(1). 

A 

This segment starts at PI(1) and terminates at the middle of the turn, marked 

A 

by the unit vector , WC(1). The logic for changing glideslope is restricted to 

A 

discontinuities at these unit middle of the turn vectors, WC(I). 

For a normal path generation sequence, arising from the first turn to 
the second turn, the existing algorithms will create three segments. The 

A 

first segment is the second half of the first turn, starting at WC(1) and ter- 

A 

minating at the outgoing unit tangent vector, PI(2). The second segment is a 

A 

great circle beginning with PI(2) and terminating at the start of the second 

A 

turn, the unit incoming tangent vector PI(3). The final, third segment is the 
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first half of the second turn about the second unit center of turn vector, 

A 

CR(2) associated with the turn radius, R^(2). 

By using the zero magnitude of the turn radius, associated with the 
glideslope change waypoint, as the index, we are able to retain all the above 
logic without change. Thus, 

If, R^(I) = 0. (8a) 

The existing algorithms will automatically set the required unit vectors 
equal to one another. 

PI(2I-1) = WC(I) = PI(2I) = CR(I) = WR(I) (8b) 

In this manner, the existing logic for providing for glideslope discontinuities at the 

A 

middle of each turn, the glideslope change wa 3 ?point, WR(I), will automatically 
be provided for. See Figure 2. 
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WR(I+2) WD 



Figure 1,— Path construction for a glideslope waypoint on 
a great circle. 
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II. Glideslope Change Waypoint Reconstruction Following Valid MLS Transition 

We assume that at the onset of a valid MLS, we are traversing a great 
circle path (straight line) containing an oncoming glideslope change waypoint. 
Utilizing the transition logic for the Continued Track, (Ref. 3), we set 

Utilizing existing RNAV equations, based on the MLS update, we may 
obtain the aircraft latitude, 6 (MLS), longitude, X (MLS), and altitude, 
h (MLS). We now set our first bulk storage wasrpoint data to be 


X (1) = 

X(MLS) 


a (1) = 

6 (MLS) 


h (1) = 

h(MLS) 

(9a) 

IC (1) = 

0 


VG (1) 

VG(RNAV) 


This defines the first wasqjoint to be 


A 

f <«MLS> ^ 


WR(1) = < 


(9b) 


L ®MLS> ^MLS> J 


This provides zero cross track error. The RNAV 

equations yield the velocity 


of the aircraft relative to the rotating Earth, in inertial coordinates, R„(t) . 

A 

The unit normal, WN(1) , is computed in the same fashion as in Reference 3 
for the Continued Track 

W(1) = jsa L ^ He (t) 

|wR(l) X Rg(t)l 

This guarantees zero track angle error. 

Since we have assumed that we are on a great circle and that the next 
waypoint is a glideslope change waypoint, the radius of next waypoint will 
be zero. 


9 



If 


R^a) = 0. 


(9d) 


then we set into bulk data 
R^(l) = R^(I) 

and 

X(2) = X(I+1) 

6(2) = X(I+1) 

IC(2) = 0 

VG(2) = VG(I+1) 


(9e) 


(9f) 


Before accepting the altitude of glideslope change wa 3 ^oint, h(I+l), we 
must first test to see if we are too close to the oncoming waypoint. This is 
necessary since we wish to avoid a double change in glideslope over a short 
time interval. 

A 

Given the existing normal, WN(1) and the latitude and longitude of the 
next waypoint, we have 


WR(2) = 



(5(2)) 

(6(2)) sin (X (2)) 
(6(2)) cos (X(2)^.. 



( 10 ) 


The desired location of WR(2) on the Continued Track great circle is given 
by Equations 5(a) and 5(b). 


WD 


1 

cos b 


^WR(2) - sin b WN(1) J 


(lla) 


where 


A A 

cosb = WR(2) * WN(1) 


sin b 



(11b) 
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To obtain the distance to go from WR(1) to WD , we have 


DIST = Rj, * (^sin”^ (I WR(1) x WD(1) [J) (12) 

If 

DIST s VERLIM 

(13a) 

h(2) = h(I+l) 

If 

DIST < VERLIM 

h(2) = h(l) + GRAD (I) * DIST (13b) 


A typical value for VERLIM is 1 kilometer, which is preset in data 
bulk storage. See Figure 3. 

The remainder of the bulk storage required for path reconstruction is 
identical to that given in Reference 3 for the Continued Path Construction. 

Finally, a change in the logic for calling in the next segment is re- 
quired if we have come to the end of the great circle segment and the next 
waypoint is a glideslope change wajqxjint or a great circle. 

If the time to go, TOGO, is triggered and R^(I) = 0. , then we update 
the segment counter by three, in order to skip over the three interior segments 
associated with the zero turn angle. 



m. Simulation Study of the Glideslope Change Waypoint Reconstruction 

A. Description of the Simulation Test Data 

This section contains computer generated plots of computer 
runs carried out using the FILCOMP program. The FILCOMP program 
was augmented to provide the algorithms described in Section I and n of this 
report to permit for the inclusion of interior waypoints at which glideslope 
changes occur along a straight line segment; Each run consists of three sets 
of plots. 

The first half-page in each series contains a plot of the aircraft ground 
track illustrating the original and reconstructed paths. The original waypoint 
data point is indicated by a point enclosed by a diamond. The reconstructed 
waypoint is marked by a point contained in a circle. The glideslope change 
waypoint is indicated by a point enclosed in a diamond (or a circle) and may be 
recognized by the fact that it occurs along a straight line segment of the ground 
track. The boundary limit of the MLS azimuth antenna are marked by a dashed 
line emanating from the azimuth antenna. The boundary limit of the elevation 
antenna is illustrated by a dashed line emanating from the antenna site to the 
right of the beginning of the runway. The initial point at the start of each tra- 
jectory is indicated by a point enclosed in a diamond, hi some instances this 
may occur too far away to be contained on the plot. 

Transition occurs immediately after all three MLS signals are received 
valid. The last to come in is usually the elevation signal and transition occurs 
when the ground track and the elevation boundary intersect. Since this series 
of simulation runs is designed to illustrate the path reconstruction for glide- 
slope change waypoints occurring on a straight line segment, the reconstructed 
glideslope change waypoint will be found adjacent to the original waypoint along 
the straight line segment when the vertical path slope discontinuity occurs. 

The cases investigated in this series include: 

(i) Glideslope change waypoint reconstruction on the first leg with a 
short distance to go. 
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(2) Glideslope change was^int reconstruction on the first leg with 
a distance to go greater than VERLIM. 

(3) Glideslope change waypoint reconstruction between two turns. 

(4) Glideslope change waypoint reconstruction on the last straight line 
segment leading to the touchdown. 

The second half of the first page in each series contains a plot of the 
altitude (meters) time history from the initial time to touchdown. 

The second page of each case consists of seven plots of pertinent data 
as a function of time. These consist of the following: 

(1) Glide path deviation in meters for both the true deviation and the 
estimated deviation for the particular navigation filter in use in the guidance 
loop for aircraft control. 

(2) Aircraft pitch angle in degrees for the true pitch, the measured 
pitch output of the IMU (used in the complementary filter) and the estimated 
pitch corrected for the estimated gjrro drift bias (used in the Kalman filter 
output). 

(3) Aircraft altitude rate, measured in meters per second, for both 
the true rate of climb and the estimated rate obtained by the particular filter 
supplying the control system equations. 

(4) , (5), (6) Errors in the estimate of the forward, lateral, and vertical 
coordinates of the aircraft measured in the flat Earth runway coordinate system 
for both the Kalman and the complementary filters, measured in meters. 

A 

(7) Error in the estimate of the forward velocity component, x. , 

Xxv 

in the runway coordinate system, measured in meters per second for both 
the complementary and the Kalman filters. 

The third and last page of each series contains eight plots of data as a 
function of time. These consist of the following; 

(1) Cross track error, measured in meters, for both the true CRTE and 
the estimated CRTE obtained by the navigation filter used to supply the guidance 
equations. 
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(2) Track angle error, converted from degrees to the time rate of 
change of cross track error by multiplying by the ground speed. Both the 
true track angle error and the estimated track angle error supplying the 
guidance equations are shown, measured in meters per second. 

(3) Aircraft roll angle, measured in degrees, for the true roll angle, 
the measured roll angle supplying the complementary filter, and the mea- 
sured roll angle corrected for the gyro drift bias (used in the Kalman filter 
estimate of the aircraft roll ai^le). 

(4) Error in the north complement of the wind, in meters per second, 
for both the Kalman and the complementary filters. 

(5) Error in the estimate of the west component of the wind, in meters 
per second, for both the Kalman and the complementary filters. 

(6) Difference between the true desired airspeed and the true airspeed. 

A second plot also shows the difference between the true ground speed and 
the true airspeed. The curves are mirror images of one another in the event 
the winds are zero, and differ in the presence of finite winds. 

A 

(7) Error in the estimate of the lateral velocity component, x , in the 
runway coordinate system, measured in meters per second for both the Kalman 
and the complementary filters. 

A 

(8) Error in the estimate of the vertical velocity component, , in 
the runway coordinate system, measured in meters per second, for both the 
Kalman and the complementary filters. 

B. Discussion of the Results 

Five cases are studied in this report. The first two illustrate 
the ability of the method to construct and to reconstruct a glideslope change 
on the initial straight line segment on the approach to the terminal area. The 
first case, Figimes 3(a), 3(b), and 3(c) illustrates the behavior of the aircraft 
in going from level flight to a 3° glideslope descent under the condition that 
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the VERLIM distance is only 60. 96 meters, (200 ft). With this small tolerance 
the original altitude is left unchanged, and the aircraft attempts to execute two 
changes in glideslope in a very short period of time, (1° sec). The effect is 
indicated in Plot 1, the glideslope deviation, in Figure 3(b). The same ex- 
treme behavior is indicated in Plot 3, of Figure 3(b), the time rate of change 
of altitude reaches a peak value of -10 m/sec (32.81 ft/sec). 

Case 2 illustrates the same approach conditions with a VERLIM of 
914.4 meters (3000 ft). Here, the distance to the approaching glideslope 
change waypoint is less than the VERLIM setting, and the altitude for the 
oncoming waypoint, at transition, is reset to maintain the initial zero glide- 
slope, so that only one charge in glide occurs. Plot 1 (Fig. 4(b)) shows a 
large reduction in glideslope deviation, and Plot 3 (Fig. 4(b)) shows a maximum 
altitude rate of 7 m/sec (21.4 ft/sec). Comparison of Case 1 and 2 illustrates 
that a VERLIM of approximately 1000 meters (3280. 8 ft) provides a comfor- 
table transition and glideslope reconstruction along a straight line segment. 

Case 3, Figures 5(a), 5(b),and 5(c) illustrate the ability to construct, 
and to reconstruct, a glideslope change waypoint along an intermediate straight 
line segment lying between two turns in the terminal area. Here, again, a 
VERLIM of 814.4 meters (2482.2 ft) is sufficient to provide a comfortable 
transition following the first valid MLS update along the intermediate straight 
line segment. Plot 3 of Figure 5(b) shows a maximum of 6 m/sec (18.3 ft/sec) 
vertical velocity required to achieve the change from zero glideslope to the 3° 
glideslope descent. 

Case 4, Figures 6(a), 6(b), and 6(c) illustrate that going to a larger value 
of VERLIM, 1524 meters (5000 ft) produces very little additional improvement 
in the response of the aircraft at transition to the required change in glideslope 
for landing. 

Case 5, Figures 7(a), 7(b), and 7(c) illustrate the capability of the method 
to design a more complicated landing approach pattern in which a change from 
3*^ descent to a zero (0.°) glideslope is required just prior to landing in order 
to avoid some obstacle. Here, the straight line segment along which the 
glideslope change occurs is the final leg immediately prior to touchdown. The 
method is seen to produce the required vertical path (Fig. 7(a)). 


15 


References 


1. Pines, S. ; Schmidt, S. F. ; and Mann, F. : Automated Landing, Rollout, 
and Turnoff Using MLS and Magnetic Cable Sensors. NASA CR-2907, 

1977. 

2. Pines, S. : Terminal Area Automatic Navigation, Guidance, and Control 
Research Using the Microwave Landing System (MLS). Part 2 - RNAV/MLS 
Transition Problems for Aircraft. NASA CR-3511 , Jan. 1982. 

3. Pines, S. ; Terminal Area Automatic Navigation, Guidance, and Control 
Research Using the Microwave Landing System (MLS). Part 3 - A Com- 
parison of Waypoint Guidance Algorithms for RNAV/MLS Transition. 

NASA CR-3512, Jan. 1982. 


16 



TABLE I 


VORTAC AND MLS STATION COORDINATES 


VORTAC STATION COORDINATES 


STATION LONGITUDE 
STATION LATITUDE 
STATION ALTITUDE 


40. 40316 
-27. 164894 
45.72 m 


MLS STATION COORDINATES 

AZIMUTH & DME LONGITUDE 40.25 

AZIMUTH & DME LATITUDE -77.025 

AZIMUTH & DME ALTITUDE 0. 

ELEVATION DISTANCE FROM RUNWAY COORDINATE FRAME ORIGIN 
XEL(l) = 1000., XEL(2) = 254.78, XEL(3) = .47 

RUNWAY HEADING = 30 ° 
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TABLE n 


INPUT DATA FOR WAYPOINT CONSTRUCTION 
CASES 1 and 2 


N = 5 


I 

X (I) 

6 (I) 

h (I) 

VdW 

IC 

Rt(I) 


Deg 

Deg 

m 

m/sec 


m 

1 

-77. 1453838 

40. 29759451 

840. 03 

74. 594 

0 

0. 

2 

-77. 1394011 

40.27189941 

840. 03 

74. 594 

0 

2286. 

3 

-77. 13148869 

40.23788553 

640.811 

69.449 

0 

1524. 

4 

-77.0584522 

40. 20574384 

290. 748 

64. 305 

0 

— 

5 

-77. 02320521 

40. 25237254 

0. 

64.305 

0 
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TABLE m 


INPUT DATA FOR WAYPOINT CONSTRUCTION 
CASES 3 and 4 


N = 5 


I 

X (I) 

8 (I) 

h (I) 

Vd(I) 

IC 

Rt(I) 


Deg 

Deg 

m 

m/sec 



1 

-77. 14538380 

40. 29759451 

490. 118 

74. 594 

0 

2286. 

2 

-77. 13148869 

40. 23798553 

490.118 

74. 594 

0 

0 

3 

-77. 10196103 

40.2249002 

490. 118 

69.449 

0 

1524. 

4 

-77. 05845224 

40. 20574384 

290. 748 

64. 305 

0 

— 

5 

-77. 02320521 

40. 25237254 

0. 

64.305 

0 

.i_ 
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TABLE IV 


INPUT DATA FOR WAYPOINT CONSTRUCTION 

CASE 5 


N = 5 


I 

X (I) 

6 (I) 

h (I) 

Vd(I) 

IC 

Rt(I) 


Deg 

Deg 

m 

m/sec 



1 

-77. 1453830 

40. 29759451 

762. 00 

74. 594 

0 

2286. 

2 

-77. 13148869 

40. 23788553 

548.64 

74. 594 

0 

1524. 

3 

-77. 05845224 

40,20574384 

195.07 

64. 305 

0 

0. 

4 

-77. 04500673 

40. 22353989 

195.07 

64.305 

0 

— 

5 

-77.02320521 

40. 25237254 

0. 

64.305 

0 

-- 
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Figure 3(b). - Case 1 
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Figure 4(c). - Case 2 . 
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